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ABSTRACT 

Context. Disclosing the structure of disks surrounding Herbig AeBe stars is important to expand our understanding of the formation 
and early evolution of stars and planets. The first astronomical units of these disks in particular, because they are hot, dense, and 
subject to intense radiation field, hold critical clues to accretion and ejection processes, as well as planet formation in environment 
different than what prevailed around our own early Sun. 

Aims. We aim at revealing the sub-AU disk structure around the 10 Myr old Herbig Be star HD 100546 and at investigating the origin 
of its near and mid-infrared excess. 

Mettiods. We used new AMBER/VLTI observations to resolve the K-band emission and to constrain the location and composition of 
the hot dust in the innermost circumstellar disk. Combining AMBER observations with photometric and MIDI/VLTI measurements 
from the litterature, we revisit the disk geometry using a passive disk model based on 3D Monte-Carlo radiative transfer (including 
full anisotropic scattering). 

Results. We propose a model that includes a tenuous inner disk made of micron-sized dust grains, a gap, and a massive optically thick 
outer disk, that successfully reproduces the interferometric data and the SED. We locate the bulk of the K-band emission at ~0.26 AU. 
Assuming that this emission originates from silicate dust grains at their sublimation temperature of I500K, we show that micron-sized 
grains are required to enable the dust to survive at such a close distance from the star. As a consequence, in our best model, more 
than 40% of the K-band flux is related to scattering, showing that the direct thermal emission of hot dust is not always sufficient to 
explain the near-infrared excess. In the massive outer disk, large grains in the mid-plane are responsible for the mm emission while 
a surface layer of small grains allows the mid and far infrared excesses to be reproduced. Such vertical structure may be an evidence 
for sedimentation. The interferometric observations are consistent with a disk model that includes a gap until ~ 13 AU from the star 
and a total dust mass of ~0.008 lunar mass (~ 6.10^^ g) inside it. These values together with the derived scale height (~ 2.5 AU) and 
temperature (~220 K) at the inner edge of the outer disk (r=13 AU), are consistent with recent CO observations. 
Conclusions. 

Key words. Stars (individual): HD 100546 - Techniques: Infrared interferometry 



1. Introduction 

■ Extrasolar planets have been discovered at a regular pace since 
] 1995. Late in 2008, five of them were imaged directly, two 
, of them clearly being associated with t he disk they formed in 
' dKalas et al.ll2008l: iLagrange et al.ll2008l) . Yet, the exact process 
[ by which they form, and the timing of their formation within 
. the disks remain unclear. There has been several interesting ad- 
' vances in the last years: direct imaging of debris disks, the dis- 
covery of fine structure in young disks likely due to gravitational 
interactions with bodies within or outside the disk, and finer 
sampling (spectroscopy) of the mid-infrared spectral energy dis- 
tribution (SED). It is now clear that several processes are simul- 
taneously acting on the disks at a given time. In the same time 

interval, it is expected that planets are actively for ming. 

In a series of multi-sca l e panchromatic stud ies dPoucet et al.l 

l2007t iPinte et all l2008bt iGradv et all l2009h . we have initi- 
ated the study of several well documented disks surround- 
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* Based on observations collected at the VLTI (ESO Paranal, Chile) 
with program 076.C-0851(B) 



ing pre-main-sequence stars. They will provide a testbed to 
study the conditions in which planets form. One of these 
stars is the late Herbig Be star, HD 100546 (B9.5Ve, 103pc; 
Ivan den Ancker et a l] (fl997h ). Its circumstellar environment was 
first resolved in scattered light by iPantin et"aD (l2000h . re- 
vealing the presence of a disk with an inclination of ~50°. 
Coronagraphic imaging showed that the star is surrounded by 
a large scale envelope (~1000 AU) and a disk extending up to 
515 AU whose brightnes s profile is asymmetric dAugereau et alj 
I2OOII: lOradv et al. 1200 lb . Its position in the HR diagram gives 
an estimated age larger than 10 Myr, supported by the presence 
of la rge fractions of crystal hne silicates that suggest an evolved 
disk dBouwman et al.ll2003h . 

The SED of the system has been intensively studied. 
Compared to other HerbigAeBe stars, its specificity is to show 
a strong mid-infrared excess, peaking at 40 jim, compared to 
a rather wea k 2-8 //m emiss i on. To account for the deficit of 
NIR excess, Bouwman et al ] d2003h suggested the presence of 
an inner cavity (empty region) within 10 AU, while fitting 
the SED with three optically thin concentric shells of different 
grain populations. Besides, support for a region of low den- 
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Table 1. Log of the observations, obtained on June 24, 2006, 
with the UT3-UT4 basehne 



Observation 


Projected 


Position 


Spectral 


Resolution 


number 


length (m) 


angle (°) 


range (fim) 


(mas) 


# 1 


59.8 


99.5 


2.06-2.46 


7.1-8.5 


#2 


60.4 


104.5 




7.0-8.4 



sity in dust grains is provided by the low circum stellar polar- 
izatio n in comparison with younger Herbig stars dClarke et al.l 
11999 *). Spectroscopic observations also agree with this conclu- 
sion. STIS observations of H2 emission re vealed the presenc e 
of a cent ral cavity extending to 13 AU (iGradv et al.l l2005l) . 
Similarlv. lAcke & van den Anckeil(l2006h concluded to the pres- 
ence of a gap distant of 6.5 AU from the star by studying the 
6300 A [OI] line. The last authors speculated that this gap (or 
secondary wall) is induced by a very low-mass stellar or plan- 
etary companion (20 Mjup). This could be the explanation for 
the pericenter asym metry observ e d in th e reflection nebulosity 
within the cavity by iGradv et"an (l2005h that could else not be 
due to chemistry effects (asumming that the density distribution 
is ax i symmetric). The recent CO observations ( Britt ain et al.l 
l2009t Ivan der Plas et al.l l2009l) reported an extended emission 
with a large inner hole inside ~11-13 AU, confirming the sta- 
tus of transitional disk given to HD 100546. 

Nevertheless, the measured NIR ex cess as well as the pres- 
ence of gas located from 0.5 to 60 AU (lAcke & van den Anckerl 
l2006t iMartin-Zaidi et al.ll2008l) . imply that the inner AUs are 
not completely empty. In the mid-infrared wavelength range, 
the disk is resolve d with mid-infrared spectro-interferometry 
dLeinert et al.l 1200 4) and nulling interferometry showed that 
the radial temperature law in the few inner AUs (at 10 fim 
and 20 jum) is n ot compatible with a continuously flared disk 
dLiu et al ] l200l . l\^nkovia et al.l (12006 ^ successfully reproduced 
the NIR SED using an optically thick disk with an inner radius 
at 0.45 AU, surrounded by an optically thin dusty enveloppe. 
All these results point toward structural differences between the 
inner and outer regions of the HD 100546 disk. 

Constraining the spatial structure and physical conditions in 
the innermost disk is therefore a crucial issue. In particular, the 
exact position and geometry of the inner rim remain poorly con- 
strained. NIR interferometric observations can achieve the spa- 
tial resolution needed to resolve this hot circumstellar material. 
The combination of angular and spectral resolution provided by 
AMBER at the ESO/VLTI allow us to direcdy resolve the cir- 
cumstellar material located in the first AU of the disk and to 
address the question of its geometry. In this paper, our goal is 
to revisit the disk geometry with a special emphasis on its in- 
ner AU. The radiative transfer models we have calculated aim at 
reproducing both the interferometric and photometric observa- 
tions. 

The paper is organized as follow: in Sect. 2, we first describe 
the observations and the data reduction procedure; in Sect. 3 and 
4, we present the interferometric results followed by the descrip- 
tion of our disk models. Finally, we discuss our results in Sect. 5. 

2. Observations and data processing 

HD 100546 was observed at the Very Larg e Telescope 
Interferometer (VLTI; iHaguenauer etal.l (I2OO8)) with the 
AMBER instrument on 24 April 2006. AMBER simultaneously 
combines three beams in the near infrared H and K bands us - 
ing spatial filtering and spectral dispersion dPetrov et al.ll2007l) . 



Table 2. Stellar and cahbrator properties 



Star 


V 


K 


H 


Sp. Type 


Diameter 












[mas] 


HD 100546 


6.7 


5.4 


5.9 


B9Vne 




HD101531 


7.3 


4.2 


4.3 


Kim 


0.6+0.1 



It therefore provides spectrally dispersed visibilities and closure 
phases. In the following, we present observations that were ob- 
tained with three 8-m unit telescopes (UTl, UT3, UT4) and 
AMBER in its low spectral resolution mode (R~30). In ad- 
dition to HD 100546, a calibrator target of known diameter 
(HD 101531) was observed. Its properties can be found in 
Table |2] The observations were achieved under good weather 
conditions (seeing of ~0.8" and coherence time of ~5ms) and 
without fringe tracker The data consist of respectively, 8 and 5 
exposures of 1000 frames each, for HD 100546 and its calibrator 
The data reduction was performed following standard proce- 
dures described in Tatulli et al. (2007), using the amdlib pack- 
age, release 2.99, and the yorick interface, provided by the 
Jean-Marie Mariotti Center (JMMC). Raw spectrally dispersed 
visibilities and closure phases are extracted for all the frames of 
each observing file. A selection process among them is necessary 
to avoid the effects of instrumental jitter and unsatisfactory light 
injection. Various selection procedures were applied, including 
selection based on the fringe signal-to-noise ratio (SNR) and on 
the fringe optical path difference (OPD; i.e. piston). Although no 
difference in absolute values was found between them, the best 
accuracy was obtained with a selection of 20% of the best expo- 
sures based on the fringe SNR. This criteria takes into account 
both the flux level and the fringe contrast (or visibility). Among 
all the measurements obtained on the three baselines, only the 
ones corresponding to the shortest baseline (UT3-UT4) in K- 
band, led to SNR above 1 .5. Consequentely, all H-band measure- 
ments, K-band data corresponding to UT1-UT3 and UT1-UT4 
baselines, and all the closure phases, were rejected. The remain- 
ing files have SNR ranging from 2.3 to 5.8. The same selection 
was applied to both the scientific star and calibrator measure- 
ments in order to guarantee an uniform data reduction. The at- 
tenuation in fringe contrast due to addition al piston with respec t 
to the zero OPD was corrected following iTatulli et al.l d2007l) . 
Finally, calibration of the AMBER-i-VLTI instrumental transfer 
function was done using measurements of HD 101531, after cor- 
recting for its diameter The errors on the calibrated spectral vis- 
ibilities include the statistical error obtained when averaging in- 
dividual frames, the error on the diameter of the calibration star, 
and the dispersion of the transfer function along the calibrator 
measurements. The last element is the one that dominates the 
error budget. After data processing, the final data set consists in 
2 measurements on UT3-UT4 baseline, that are spectrally dis- 
persed in 13 K-band spectral channels. Therefore, in total, we 
obtained 26 visibilities in the [2.06-2.46] jum wavelength range. 
A summary can be found in Table [1] 

3. Results 

3.1. Visibilities 

Fig. [U left, shows the measured visibilities after data process- 
ing. They are plotted against the spatial frequency, B/A, with B 
the projected baseline length and A, the wavelength of the ob- 
servations. The spatial frequency is the inverse of the resolution 
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Fig. 1. Left: measured visibilities plotted against spatial frequency. The two different measurements (#1, #2) are shown with distinct 
symbols. Right: the visibilities (black crosses) are shown against wavelength and are compared to the predictions of a ring of uniform 
brightness. The red full and dashed lines give these predictions for the two measurements {i.e. at two distinct projected baselines). 



achieved by the interferometer and takes into account the wave- 
length, the projected baseline, and consequentely, the direction 
on the sky. It can be expressed in unity [m/yum] or in cycles per 
milli-arcsecond. It can be immediately seen that the circumstel- 
lar matter around HD 100546 is resolved at the level of a few 
milli-arcseconds, because the visibility is well below 1. 

To derive basic characteristics of the NIR emission, one can 
consider a simple geometrical model of a ring of uniform bright- 
ness that roughly accounts for hot dust located in the inner rim 
of a dusty disk. The model visibilities Vmod can be written : 

_ K + Vex * F,JF, 

with V*=l (unresolved central star) and Fex/F,, the ratios of the 
NIR excess to the stellar flux. This ratio, ~2.2, is estimated at 
each AMBER wavelength, using a Kurucz model for the central 
star emission at T= 10500 K. With this model, the extension of 
the K-band emission is found to be ~5.4 mas (i.e. a radius of 
~0.28 AU) at d=103 pc. The observed chromaticity of the visi- 
bilities can be compared to the predictions of this model. Fig.[T] 
right, shows the spectral visibilities against wavelength. The vis- 
ibility curve is almost flat with wavelength, while the uniform 
brightness ring predicts a slightly steeper slope. 

HD 100456 was also observed by Leinert e t'aD (|2004 in the 
mid-infrared with the MIDI instrument at VLTI. These observa- 
tions led to a single-dish image at 8.7 fim and to visibility mea- 
surements (at 74 m baseline length). In Sec. 4, we use this image 
and the corresponding 8.7 fim visibility. 

3.2. Photometry 

We retrieve from the litterature photometric measurements from 
the UV to the millimetric wavelengths t o produce a complete 
SED (2MASS, IRA S, ISO ; lArdila et al.1 (l2007h : iHenning et all 
(ll994l) : lGradv et all (l2001h ). Moreover, because young sta rs are 
highly variable in the NIR wavelengths (" Sifl^io et al. 1'2008'). and 
since modelling of the observations relies on these photometric 
values, it is useful to provide an estimate of the NIR emission 
directly from the interferometric measurements. We were able 



to extract the K band photo metry from the thre e AMBER photo- 
metric channels. Following I Acke et al.l (12008'). for each file, we 
summed the flux in each channel over all frame and calibrated 
this value with the total flux in the corresponding channel for the 
calibrator We applied a photospheric model for the calibrator 
star in order to convert it in flux units and obtain 5.33+0.08 as 
the K band magnitude, consistent with the 2MASS values. 



4. Modelling 

Since it ha s been long estab lished that HD 100546 is surrounded 
by a disk (lAugereau et al.l 12001.) . in this section, we use disk 
models to describe its circumstellar environment. It is not pos- 
sible to unambiguously estimate the disk inclination and posi- 
tion angle from a single-baseline interferometric measurement, 
therefore, we adopt the previously pu blished values of 42° and 
145°, respectively (lArdila et al.ll2007l) . Besides, we consider a 
total disk mass of 5 x10''* M©, inferred by 1.3 mm imaging 
(IHenning etal.lll994 . 

4.1. The radiative transfer code MCFOST 

We used the Monte Carlo-based 3-D radiat ive transfer code , 
MCFOST, to produce disk images and SED (iPinte et al.ll2006l) . 
The temperature distribution is calculated using the immedi- 
ate re-emission concept of Biorkman & Wood (2001) but with 
a co ntinuous d eposition of energy to estimate the mean inten- 
sity (lLucvlll999 ). In this paper, we used a cylindrical grid, with 
an adaptive mesh at the inner edge (based on the opacity gradi- 
ent) so as to properly sample the inner radius of the disk. The 
radiation field and temperature structure estimated by the Monte 
Carlo runs are used to produce SEDs by calculating the formal 
solution of the radiative transfer equation along rays (ray-tracing 
method), as well as images. Visibilities are computed by taking 
the Fourier transform of the latter, at the specific spatial frequen- 
cies sampled by the observations. The code has been extensively 
tested and its performan ces compared to other radiative transfer 
codes (iPinte et al.ll2009l) . 
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Astronomical Units Baseline (m) Wavelength (microns) 



Fig. 2. Left: K-band model image of HD100546. The full and dashed lines indicate the position angles of the observations (PA = 
99.5° and PA - 104.5°, respectively). Middle: model visibilities, as a function of baseline length. The broad band K-band AMBER 
visibilities are overplotted. Right: spectrally dispersed AMBER visibilities versus wavelength, for both measurements (upper and 
lower pannels respectively). The spectral profiles of the model visibilities, as derived from our best model, are added. 



We assume that the disk is constituted of silicate dust grains 
only (Pollack et al. 1994), and is in hydrostatic equilibrium. Its 
structure is parametrized by power laws, so that the surface den- 
sity 2, and the scale height H can be expressed as varying with 
radius, r, following : 

£(r) = Z,ooAu(y^)' 

H(r) = H,„„Au(^^f 

j6 is the flaring exponent. Siooau and Hiooau are the values of 
the surface density and the height scale at a radius of 100 AU. 
The dust grain sizes, a, are distributed following a power law 
a~^'^ (ISM), with a varying between the minimum and maximum 
sizes, a^in and a^ax respectively. In addition, we fix the total dust 
mass (Mdust) of the disk, that is delimited by its inner and outer 
radii (Rin and Rout)- 

The disk of HD 100546 is well known and i ts geometry 
and content have been estiniated by several autho rs (iGradv et al.l 
l200lHBouwman et al.ll200llBrittain et alJl2009( ). We will there- 
fore not provide a detailed re-analysis of all the disk structural 
parameters. To the contrary, using in spirit the results of previous 
studies, we will construct a disk model that also includes a large 
gap. However, with the help of interferometric measurements, 
and complete radiative transfer including scattering, we will re- 
visit the morphology and mass content of the inner disk, as well 
as the position of the gap, in a global picture. 

4.2. The disk model 

The SED constrains the mass-temperature law for the dust. 
The near-infrared emission in excess of the photosphere is 
due to hot material at ~1500 K, and must result from a small 
amount of dust in an optically thin inner disk. The large mid 
and far-infrared excesses are due to dust at a temperature of 
~200 K, while large and cold grains in the mid-plane of the disk 
must be responsible for the emission at longer wavelengths. 
The SED gives integrated flux informations but provides little 
morphological informations on the dust spatial distribution. 
In this section, we propose a disk model that successfully 
accounts for the emission at all wavelengths, as well as for the 



AMBER and MIDI visibihties (see Fig. |2] |4]and|5]l. However, 
our model does not attempt to match mineralogical features. As 
pointed out by iBouwman et al.i (i20()3l) . continuous flared disk 
models can not reproduced the SED as they produce either too 
much or too little flux in the near and mid-infrared wavelengths, 
respectively. The model includes an optically thin inner disk, a 
gap and an optically thick outer disk. The model parameters can 
be found in Table |3] including for each zone, the mass (Mdust), 
the inner and outer radii (Ri,,, Rout respectively), the scale height 
at 100 AU (Hiooau), the flaring and surface density exponents 
(J5, q, respectively) as well as the minimum and maximum grain 
sizes (amin, ^m^x, respectively). 

The inner disk: The AMBER/VLTI measurements can 
locate the bulk of the emission in the K-band at ~0. 25-0. 30 AU. 
This puts a strong constraint on the location of the disk inner 
edge and consequently on the silicate dust grains sizes. The 
tenuous inner disk has a total dust mass of 3.10"'° Mq and starts 
at a radius of ~0.26 AU, a distance inferred from sublimation 
of micron-sized silicate grains (specifically, with sizes ranging 
from 0.1 to 5 fim) at ~1500 K. Its outer radius is not well 
constrained by the AMBER observations because it is fully 
resolved, and we fix it at ~4 AU. Here, all the model calcu- 
lations are done assuming that the dust grains are at the same 
equilibrium temperature with the stellar radiation field which 
may not be exactly the case in the tenuous optically thin inner 
regions. Although a more detailed calculation, that requires 
determination of the dust-to-gas ratio is beyond the goal of this 
paper, we stress that there may not be a specific inner radius but 
a range of inner radii, depending on the grain sizes. 

Tlie gap: After computing numbers of models, we confirm 
the need for a very low density region [e.g. a gap) to reproduce 
the lack of emission in the 3-10 fim wavelength range. The 
MIDI/VLTI measurements provide a complementary informa- 
tion on the spatial distribution of the dust in the disk. At the 
achieved angular resolution at ~8.7 jum (24 mas i.e. 2.4 AU), 
the innermost disk is not resolved but the outer disk is fully 
resolved. The 8.7 yum observations therefore bring a strong 
constraint on the flux ratio, between these two elements of 
the disk. In the model, the gap starts at 4 AU until 13 AU, a 
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~ 350 AU 

Fig. 3. Schematic view of our best disk model, that includes different zones: an inner disk with micron-sized silicate grains, a gap, 
a massive outer disk with small grains in an upper layer. 




X (\im) X (ixm) 

Fig. 4. Left: measured SED (circles) for HD 100546, compared to the predictions of our best model (black full line). The individual 
contributions of each disk zone to the SED are shown (inner disk, red dashed curve; surface layer, blue dotted curve; outer disk, 
green dash-dotted line). A similar fit to the SED can be obtained with a diflFerent disk model (magenta line). Right: zoom in the near 
and mid-infrared wavelengths. 



radius consistent \y ith spectroscopic results (iGradv et al.ll2005t 
iBrittain et al.]|2009l) . 

The outer disk and its surface layers: the model massive 
outer disk forms at 13 AU and extends up to 350 AU. It contains 
most of the disk mass and holds large dust grains (from 1 yum to 
10 mm). To account for the mid and far-infrared excesses, we 
include a population of small grains {i.e. 0.05 to 1 fim) at the 
surface of the outer disk, at distances of a few tens of AUs (here, 
from 13 to 50 AU; see Table |3]l. The latter are not precisely 
constrained by the observations, as the only requirement is to 
have a specific amount of small grains emitting at a few tens of 
AUs. 

The entire set of parameters used in this model is given in 
Table [3] and a schematic view of the disk model is provided 
in Fig. [3] In Fig. |4] we report the contributions from each part 
of the disk to the total flux. The inner disk (red dashed line) 
produces enough emission to fit the NIR excess while the small 



grains in the surface layers (blue dotted line) produce the major 
contribution to the large mid and far infrared excess. Finally, 
the large grains in the outer disk are responsible for the flux at 
long wavelength (green dash-dotted line). The predictions of the 
entire model are in good agreement with the observations (full 
black line). 

For this model, images and visibilities were computed at the 
AMBER wavelengths. Figure |2] (middle) shows the visibility 
profile as function of the projected baseline length. This model 
succesfully reproduces the interferometric observations in addi- 
tion to the SED: the spatial extent of its K-band emission and its 
inner structure are in good agreement with the data. The chro- 
maticity in the K-band is also well reproduced with the radial 
temperature profile of the disk model (Fig.|2] right). In the mid- 
infrared, we focus our modelling on the continuum emission at 
8.7 yum and do not attempt to reproduce any spectral feature. 
Figure |5] (left) gives the model image at 8.7 jjm that clearly 
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Fig. 5. Left: 8.7 /im continuum image of HD 100546 predicte d by our best m odel and oriented along a 145° position angle. Middle: 
8.7 jum continuum image obtained with a sing le UT at VLT dLeinert et al.ir2004,.) . For both images, the contour levels are linearly 
spaced from 100% to 25% as the fourth root of the image brightness. R ight: visibility curve calculated for our best model (solid 
line) at 8.7 jjm. The MIDI measurement obtained bv lLeinert et alj (|2004|) is shown with the error bar. 



Table 3. Parameters of our best model. 





inner disk 


surface layer 


outer disk 


Mdust 


[Mo] 


3e-10 


6e-5 


5e-4 


Rin 


[AU] 


0.26 


13 


13 


Rout 


[AU] 


4 


50 


350 


HloOAU 


[AU] 


6 


12 


12 






1 


0.5 


1.125 


q 




-1 


-1 


-1 




[jum] 


0.1 


0.05 


1 


^max 


[/jm] 


5 


1 


10000 



Table 4. Contributions from the different types of emission to 
the total flux. 



Model with 


Fdir.-st/Ftot 


Fsc.st/Ftol Fdir.lh/Fiot 


F.sc.ih/Ftoi 


Ri„=0.26 AU 


19.6 


4.1 38.0 


38.3 


Ri„=0.40 AU 


16.7 


1.8 64.4 


17.1 



shows an e longated disk struct ure, in agreement with the one 
obtained bv lLeinert et al.l (l2004l) (Fig.|5] middle). From a simple 
2-D fit, we obtain ~21 AU and ~14 AU as the disk major and mi- 
nor axis, respectively, in good agreement with the values derived 
by Leinert et al. 2004 (29 +4 AU and 19 +9 AU respectively) and 
compatib l e with the results from nulling interferometry (24 AU; 
iLiu et al ] (l2003h 1.Figurel5l right, shows the model visibilities at 
8.7 yum, compared to the MIDI measurement. At long baselines, 
the outer disk is entirely resolved and one has access to the unre- 
solved flux contribution of the inner disk to the total 8.7 //m-flux. 
At a 74 m-baseline, our model can well account for this flux ra- 
tio. Our model predicts that, at short baselines (<10 m), the outer 
disk is resolved as the steep decrease of visibility attests, while 
for intermediate baselines, the sharp edge of the gap produces 
a clear signature. Mid-infrared interferometry therefore appears 
well suited to probe the presence of a gap, as the strong gap sig- 
nature can easily be detected using moderate VLTI baselines. 

5. Discussion 

In this paper, we propose a model of the disk surrounding 
HD 100546 that successfully reproduces the SED, AMBER 



and MIDI interferometric observations. We note that equally 
good fits to the SED can be obtained with other models, similar 
in spirit to the one we propose (i.e. with three distinct disk 
regions) but wi th different spatial distribution. In particular, 
[Bouwman et al,l |2003) suggested a solution to reproduce the 
SED, that includes a population of small grains distributed from 

0. 3 to 9.8 AU, a more massive part with small grains from 9.8 to 
43 AU, and a population of large dust grains from 28 to 380 AU, 
with a grain size distribution n(a) oc a"^. Different species are 
used to study the mineralogy of the dust, including refractory 
grains, such as iron. Using the same grain sizes distribution, but 
for silicate only, we could produce a good fit to the SED (see 
Fig. m left, magenta curve). To account for the NIR emission, 
the small grains are nevertheless required to be hotter than 
the tabulated evaporation temperature for silicate (i.e. 2100 K 
instead of 1500 K). Besides, they induce a slightly different 
location for the bulk of the K-band emission (i.e. 0.4 AU) which 
does not match the AMBER visibilities. We therefore confirm 
here that spatially resolved informations are needed to speculate 
on the disk structure and properties, as fitting the SED alone is 
highly ambiguous. 

The importance of scattering: An interesting side effect 
is the amount of scattered light induced in different models. At 
each wavelength, the total emission (Ftot) is constituted of 4 
distinct contributions: the stellar light that directly reaches the 
observer (Fdir,st); the stellar light that is scattered by the disk 
(Fsc.st); the direct thermal emission from dust grains (Fdir,th; 

1. e. absorbed and then re-emitted photons) and finally, the 
scattered thermal light (Fsc,th; i-e. absorbed, re-emitted and 
then scattered photons). Micron-sized silicate grains located 
at 0.26 AU are good scatterers, since their sizes are of the 
order of the AMBER wavelengths. In our model, the scattered 
thermal light contributes to -38% of the total K-band emission, 
about as much as the direct thermal flux. On the contrary, the 
model with small grains in the inner disk and an inner radius at 
0.4 AU, predicts a dominant contribution from direct thermal 
light (~64%). Table |4] reports the individual flux contributions 
to the total K-band emission (Ftot). The index 'dir' or 'sc' 
indicates whether it is direct or scattered light, while the index 
'st' or 'th' designate stellar and thermal light, respectively. 
Assuming that the disk is made of silicate grains only, the fit 
to the AMBER data constrains the grain size and consequently 
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the amount of scattered light. This rules out solutions where 
only direct thermal emission from the dust is assumed. This 
means that scattered light can have a major effect in the overall 
energy budget in the NIR, and a strong impact on the inner 
disk parameters inferred from modelling. This was already 
showed and taken into account in modelling obse rvations of 
solar-mass young stars, as previously shown by iPinte et aP 
(2008a), but this is the first time that such case is made for a 
Herbig Be star This conclusion in fact contrasts with the usual 
assumption made to interpret NIR interferometric results in 
Herbig AeBe stars, i.e. that the NIR excess emission is due 
only to direct thermal emission from hot dust located at the 
inner edge of a disk ( Dullcmoiid el a l. 200 1 ). In fac t, scattering 
or a dditional emitting compo nents (such as gas dlsella et alj 
12008 ) or refractory dust grains (i Benistv et al.ll2009l) ) seem to be 
needed to reproduce recent spectro-interferometric observations 
of Herbig AeBe stars. Whether this is relevant for all the objects 
of this class should be validated on simpler cases. 

Evidence for sedimentation? The need for small grains at 
the surface layers and for large grains in the bulk of the disk can 
be interpreted as resulting from grain growth and dust settling. 
A toy model was computed with dust grain sizes ranging from 
0.05 to 3000 fim, a standard flaring (j6= 1.125), but considering a 
larger scale height (36 AU at 100 AU) and a strong dust settling. 
The predicted emission similarly fits the SED. Considering the 
age of the system, the disk is very likely evolved and vertical 
stratification could well be happening. 

A gap and a tenuous inner disk: Modelling of the SED re- 
quires the presence of a gap. In fact, a continuous disk model 
can fit the AMBER data and the flux at 2.2 fj.m, but it produces 
too much flux at 8.7 jum. Similarly, a continuous disk model that 
would fit the MIDI data and the 8.7 fj.m flux, would in turn pre- 
dict too little flux in the NIR. This confirms the need for a gap, 
or at least, for a strong discontinuity in the radial density pro- 
file that would produce the needed amounts of flux at specific 
locations. 

If the need for this density discontinuity is clear, its size is 
nevertheless very model-dependent. In fact, we can not deter- 
mine the position of the outer radius of the inner disk with preci- 
sion, because it is fully resolved by AMBER. The MIDI obser- 
vations constrain the flux ratio between the inner and outer disks: 
when our inner disk extends up to ~4 AU, the outer disk starts 
at ~13 AU. A change in the inner disk geometry would affect 
this estimate. In our model, it is constrained by the temperature, 
hence position, of the emitting particles. Considering the typi- 
cal grain size distributions given in Table [3] our estimate is in 
good agreement with the values inferr ed from previous studies 
ferittain et alJl2009HGradv et alJl2005h . 

iBouwman et alJ (l2003h required the disk to be more puffed 
up at 10 AU than prediced by a standard flared disk at the same 
radius to capture and re-emit enough light. They use 3.5 AU as 
scale height for dust at ~200 K which is much more than ex- 
pected for gas at the same temperature (i.e. 1 AU). 'Britta in et al] 
(2009) found that the CO disk extends from 13+6 AU and that at 
the inner edge of the disk, the gas is much hotter (~1400 K) than 
the dust at the same radius. As the authors suggest, this can be 
explained if there is a low density region in the upper disk atmo- 
sphere, that results in a scale height of ~3 AU. In our model, the 
scale height of the outer disk at its inner edge is ~2.5 AU, due to 
the small grains i n the surface layers, a value in agreement with 
the predictions of lBrittain et alJ (l2009l) . 



Moreover. lBouwman et"a n (l2003h predicted a total dust mass 
of 10^^ g inside the gap. This led 'B rittain et al.1 (l2009l) to de- 
duce from their estimate of the mass of gas in the inner disk 
(<3xlO^^ g) that the gas-to-dust ratio is ~3. Our model re- 
quires a total dust mass of 3x10""^ Mq (6x10^-^ g) i.e. a gas- 
to-dust ratio of ~50. We therefore conclude that these esti- 
mates are highly model-dependent and should be considered 
with care, pending a direct measurement of the dust content of 
the inner disk. The small accretion rate, inferred from UV spec- 
troscopy dDeleuil et al.ll2004t iGradv et al]l2005h - a few times 
lO"^M0.yr~'- is enough to empty the cavity in a few decades, 
yet we are observing dust in the inner disk at an age >10 Myr 
Either we are observing HD 100546 at a very special moment 
in its evolution or the inner disk is being replenished. We note 
that a putative large orbiting body, a Jupiter-like planet located 
at about 10 AU, has the capacity to carve the disk and lead to 
such a surface density discontinuity at ~13 AU, while allowing 
some mate rial to flow in contin uously across the gap all the way 
to the star (IVarniere et alJl2006h . 

This transitional disk will be a primary target for ALMA 
that will enable direct determination of the gap extent, by ac- 
curately probing the gas surface density profile. Before then, ad- 
ditional VLTI data with three telescopes, allowing phase closure 
detections, will add new constraints such as the asymmetry of 
the inner disk rim, the co-planarity of the inner and outer disks 
(Gradv et al. 200 9), and whether t he star is off-center with re- 
spect to the gap (iGradv et alJl2005h . 
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